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We consider the helical reduction of the wave equation with an arbitrary source 
on (n+ 1) -dimensional Minkowski space, n > 2. The reduced equation is of mixed 
elliptic-hyperbolic type on R n . We obtain a uniqueness theorem for solutions on 
a domain consisting of an n-dimensional ball B centered on the reduction of the 
axis of helical symmetry and satisfying ingoing or outgoing Sommerfeld conditions 
on dB m S n . Non-linear generalizations of such boundary value problems (with 
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1. Introduction 

Recent approaches to the quasi-stationary approximation to the intermediate phase of 
binary inspiral in general relativity have led to the consideration of reductions of the Ein- 
stein equations by a helical Killing vector field (see [1,2,3] and references therein). To date, 
model problems have been analyzed consisting of helical reductions of linear and non-linear 
wave equations in (3+l)-dimensional Minkowski spacetime with various sources using Som- 
merfeld conditions on a spherical boundary. These helically-reduced equations have the 
challenging feature of being of mixed elliptic-hyperbolic type on their 3-dimensional do- 
main. More precisely, they are elliptic in an inner cylindrical region surrounding the sources 
and hyperbolic outside this cylindrical region. There appear to be no general theorems 
to handle existence and uniqueness of solutions to partial differential equations of mixed 
type. Results tend to be specific to individual equations or limited classes of equations, 
and even then the equations which have been most studied are defined in 2 dimensions 
[4] . From the investigations of [1] it appears that the boundary value problem arising from 
helical reduction of (linear and non-linear) wave equations using Sommerfeld conditions 
on an exterior boundary is well-posed. Solutions have been constructed and appear to be 
unique. This is somewhat remarkable since the boundary intersects both the hyperbolic 
and elliptic domains. In particular, one might not expect a single (Sommerfeld) condition 
on a closed boundary to enforce uniqueness of solutions [1]. 

Some light was shed on this issue by the work of [5] where the helical reduction of 
the (2 + l)-dimensional wave equation was shown to define a symmetric-positive system 
on an annular region in R 2 such that the Sommerfeld boundary value problem was well- 
posed — solutions exist and, in particular, are unique. Unfortunately, it is not known how 
to generalize these results (i.e., symmetric positivity of the reduced equation) to higher 
dimensions. Moreover, the helical reduction of the (2 + l)-dimensional wave equation 
leads to a boundary value problem on a two-dimensional region with an outer circular 
boundary which need never intersect the circle of degeneracy of the symbol of the reduced 
partial differential equation. In higher dimensions, the spherical outer boundary necessarily 
intersects the "light cylinder" where the symbol is degenerate so the boundary conditions 
must be imposed both in the elliptic and in the hyperbolic regions.* This makes the 
problem qualitatively different in the physical (3 + 1) spacetime dimensions (and in higher 
dimensions). 

Thus it is of interest both from mathematical physics and gravitational physics view- 
points to better understand the nature of boundary value problems arising from helical 
reduction of wave equations. Here we shall provide a uniqueness theorem for the heli- 

* Unless, of course, the boundary is completely contained in the elliptic region, which is not 
of physical interest and which, in any case, leads to a standard elliptic boundary value 
problem. 



2 



cal reduction of the (n + 1) -dimensional wave equation with arbitrary sources and with 
Sommerfeld boundary conditions. The proof is remarkably elementary and employs an 
approach used by Protter to study a generalization of the Tricomi problem [6] . 

2. The helically-reduced wave equation 

We will be considering the helical reduction of the wave equation with an arbitrary 
source on (n + l)-dimensional Minkowski space, with n > 2. The spacetime manifold is 
N = R n+1 with metric 

77 = —dt <S> dt + dx <S> dx + dy <E> dy + dijdz 1 <g> dz J , (2-1) 

where Latin indices i,j = 1, 2, . . . , n — 2. The wave equation for ty:N — > R with a 
prescribed source F: N — > R is 

+ * xx + Vyy + fiVij = F (2.2) 

Note we use the notation where subscripts on a function indicate partial derivatives. The 
helical reduction is accomplished by assuming the source and solutions are invariant with 
respect to the isometry group (G) generated by 

K = dt + fl(xd y — yd x ) , O = const., (2.3) 

which is equivalent to 

C K F = C K ^ = 0. (2.4) 
In cylindrical coordinates (£, p, 0, z l ), the metric and Killing vector field are 

77 = -dt ® dt + dp ® dp + p 2 d(j) ®d(j) + dijdz 1 <g> dz\ (2.5) 

K = d t + O^, (2.6) 

the wave equation is 

+ ^(ptfp) + + = F, (2.7) 

and the invariance condition (2.4) is 

*t = -H^, ^ = -fii^, (2.8) 

Introducing (p = <f> — fit, (2.8) means there exists functions u and / such that 

p, 0, z*) = U (p, v?, A F(t, p, <f>, z l ) = f(p, <p, z l ). (2.9) 
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We then get the reduced equation defining helically-invariant solutions to (2.2) 

- p d p {pu p ) + ^-u w + = /, (2.10) 

where 

x(p) = i - nV- (2-n) 

The locus of points where x(p) = is the "light cylinder". Inside the light cylinder 
(p < jj) eq. (2.10) is elliptic and outside the light cylinder (p > ^) eq. (2.10) is hyperbolic. 

A useful geometric interpretation of this reduction is as follows. The set of orbits of 
the group generated by K defines a manifold M = N/G R n . The functions (p, <p, z l ) are 
G-invariant and define cylindrical coordinates on M. In these coordinates the projection 
7r: N — > M is simply 

n(t,x,y,z l ) = (p,<p,z l ), (2.12) 

and satisfies n*K = 0. The G-invariant functions F and \l/ on N correspond to functions 
/ and «onM, respectively, via 

F = ic*f, * = n*u. (2.13) 
The inverse metric on iV is given by 

= -d t ® d t + dp ® d p + + 8 ij di dj . (2.14) 

Being G-invariant, rft projects to a tensor field q on M. Using (2.12), 

q = Tr^tt = dp ® dp + ^r-d v ® 8^ + fidi ® dj. (2.15) 

This tensor field is well-defined everywhere on M, but it does not determine a metric on 
M because a has no inverse on the light cylinder. While the metric on iV does not induce 
a metric on M, the metric volume form e on N does define a volume form v on M as 
follows. Define 

u = K->e, (2.16) 

which satisfies 

L K u = 0, K^u = 0. (2.17) 
Consequently, to is the pull-back by ix of a volume form v on M. It is easy to check that 

z/ = pdp A dtp Adz 1 A . . . A cfe n ~ 2 . (2.18) 



The volume form v defines a scalar density of weight-1, a = p, on M. 
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We will use Greek indices to label tensor fields on M. Introduce a torsion- free derivative 
operator V a . The reduced equation (2.10) is equivalent to 

-V a (?q afi Vpu) = f. (2.19) 

To see this, we first note that, because of the density weights, (2.19) is in fact independent 
of the choice of torsion- free derivative V«. Using the cylindrical coordinate derviative 
operator, V a = d a , in (2.19) we obtain (2.10). For what follows we re-write (2.19) as 

V a {h a N p u) = d a (h a %) = /, (2.20) 

where / = af is a scalar density of weight one and h a @ = aq a P is a tensor density of 
weight- 1 given by 

hPP = p, W=pb l \ = - x =--n 2 p. (2.21) 

p p 



3. Energy integral 

The key ingredient in our uniqueness theorem is the following generalized energy inte- 
gral. Fix a domain B C M and define 



E[u}= / \(au + b^ Ul )d a (h a ^up) 



D 



(3.1) 



where a and b a d a are a function and vector field to be specified later. The integrand 
involving a can be written as 



aud a (h a P U(3 ) = ^da(h a ^ a/3 )u 2 - ah a ^u a up + d, 
The integrand involving K 1 can be written as 

(3.3) 

Again, while these expressions use the coordinate derivative, they are in fact independent 
of the choice of torsion-free derivative operator. The divergences integrate to the boundary 
and we have 



(3.2) 



E[u] = J S^d a (h a0 a p )u 2 - ah a(3 u a u p + ^(h a ^)u a u - b^h^UjU 

+ [ natiau + b^u^h^up-lh^apu 2 -h a h^u lU 
JdB I 2 2 



(3.4) 



If there were a metric on B, n a could be defined in terms of the unit normal to the 
boundary and the metric-induced volume element of the boundary. Without a metric 
n a is still defined, of course, but its definition is necessarily more involved. We give the 
definition in the Appendix. 



4. Uniqueness theorem 

We are now ready to formulate the boundary value problem of interest. We consider 
solutions to the equation (2.20) on a ball of radius R: 

B = {(p,(p,z*)\0 < p 2 + S lJ z t z J < R 2 }. (4.1) 

The boundary dB is the sphere S 71 ^ 1 of radius R. Using (A. 7), we have in spherical 
coordinates (r, 0\, . . . , 9 n -\) on B : 

n a dx a = dr. (4.2) 

We impose Sommerfeld conditions on dB. Taking account of the helical reduction they 
are of the form 

— (pu p + z l Ui) ± Qd(pU = r, on dB (4.3) 
R 

where r: S n ~ l — ► R is some specified function. 

We remark: (i) if the boundary is chosen such that R> ^ the boundary passes through 
both the elliptic and hyperbolic domains; (ii) r and / cannot be specified independently; 
the integral of (2.20) over B implies 



f f = f or. (4.4) 

Jb job 



Since (2.20) and the boundary conditions (4.3) only involve derivatives of u, solutions 
to these equations can only be unique up to an additive constant. In fact, this is the only 
freedom in the solution. Our main result is the following. 

Theorem. Given O, f:B^ R, and r: dB — > R, any two solutions to (2.20) on B with 
boundary conditions (4.3) differ at most by a constant. 



Proof: 



Consider the difference of two solutions, u = u\ — u% u satisfies (2.20) and (4.3) with 
/ = and r = 0, respectively. Consequently, E[u] =0 for any choices of the function a 
and vector field b = b a d a . We choose these as 

2 



-1, 



1 -n 



pd p + z l d, ± Rm 



(4.5) 
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Note that 

b a u a = 0, on OB. (4.6) 
A straightforward computation, using (4.3) with r = in the boundary integral, then gives 



= 



B I \ n - 1 



(u 2 p + S ij u iUj ) + (-j + O 2 )^ 



+ 



/■ / 1 \ B r , (47) 

y ( ^1 J {(^ + P 25%J ) u i u j + (1 + ^ 2 SijZ % z 3 )ul ± 2Bnz l u iUip j . 

The volume integrand (in the first integral) is manifestly non- negative for n > 2. We now 
show that the boundary integrand (in the second integral) is also non-negative. 

We first note that the boundary integrand is invariant under orthogonal transforma- 
tions of the z l . Thus, given any point (p,(p,z l ), we can rotate the z l axes such that 
z 1 = (z, 0, 0, . . . , 0), where z 2 = 5ijZ l zK The boundary integrand at the given point is then 

^~Tl) {^"l + P 2 $ ij UiUj + (1 + ftV)u£ ± 2i2fizuiu v } ^ ^ 



> 
> 0. 



Because both integrands in (4.7) are non-negative they must each vanish. From the 
volume integrand it follows immediately that 

u a = 0. (4.9) 
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Appendix: The divergence theorem without a metric 

Consider an n-dimensional orientable manifold M, a torsion-free derivative operator 
V a on M, and a vector density of weight one V a . Given B C M, Stokes theorem implies 
an identity of the form 

/ V a U a = / n a V a . (A.l) 
JB JdB 
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Normally this divergence theorem is proved using a metric on M. However this is not 
necessary. Here we shall give a version of the divergence theorem and, in particular, give 
a formula for n a without using a metric. 

The manifold M, being orientable, comes equipped with a nowhere vanishing n-form 
density of weight minus 1, denoted by r} a - l --a n i an d a totally antisymmetric contravariant 
tensor density of weight one, fj a i"' a n^ guch that 

V n V^ n =nl6%\..6« : \ (A.2) 

Both rj av ..a n and are constant for any choice of V^. 

The boundary dB is an oriented submanifold in M embedded by i: S — > M, i.e., 
dB = i(S). S is equipped with an (n — l)-form density of weight minus one, £ai...a n _D an d 
a skew, contravariant rank (n — 1) tensor density of weight one, °n-i ) satisfying 

i ai - an -^ bl ... bn -i = (» - W [ £ ■ • • ( A - 3 ) 

(In this Appendix only we use Latin indices to denote tensors on S.) 
To apply Stokes theorem we define an (n — l)-form 

Wai...a n _i=^Vi-a n -r ( A - 4 ) 
We then have (using differential form notation) 

/ V a V a = [ dw 
JB JB 

u (A.5) 



(A.6) 



dB 

= / 5 (^l)! ri '" fl " 1(fw)ai - a "- 1 - 
Now, at points of dB we can write 

Thus we have 

An alternative approach to the integral over B in (A.l) is to note that it is independent 
of the choice of V«. If we fix a Riemannian metric g a p on M, and use the metric compatible 
derivative operator, we have available the more traditional form of the divergence theorem: 

/ V a V a = [ ^h a W a , (A.8) 

JB JOB 



where 

W a = —V a , (A.9) 

Vg 

n a is the outwardly oriented unit normal to dB, and ypy is the induced volume element on 
dB. The result (A. 8) is, of course, equivalent to the manifestly metric independent result 
(A. 5) above, as can be verified by using the identity 

H = (^I)!^^^)" 1 -""" 1 ^-^^- (A.10) 
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